Hot gets hotter through long-range contact with a cold bath 
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We investigate the dynamics of a small long-range interacting system, in contact with a large long- 
range thermal bath. Our analysis reveals the existence of striking anomalies in the energy flux 
between the bath and the system. In particular, we find that the evolution of the system is not 
influenced by the kinetic temperature of the bath, as opposed to what happens for short-range 
collisional systems. As a consequence, the system may get hotter also when its initial temperature 
is larger than the bath temperature. This observation is explained quantitatively in the framework 
of the collisionless Vlasov description of dynamics and shown to be universally valid whenever the 
Vlasov picture applies, from cosmology to plasma physics. 



PACS numbers: 05.20.-y,05.70.Ln,52.65.Ff 

In the recent past, several theoretical and experimental 
studies have been devoted to exploring dynamical and 
thermodynamic properties of long-range interacting sys- 
tems (LRIS's) pp. In such systems, energy is not additive. 
This fact, together with a possible break of ergodicity, 
is at the origin of a large gallery of peculiar thermody- 
namic behaviors: the specific heat of LRIS's can be nega- 
tive in the microcanonical ensemble [2] and temperature 
jumps may appear at microcanonical first-order phase 
transitions. These systems also display remarkable non- 
equilibrium dynamical features. For example, it is well 
known that under particular conditions isolated LRIS's 
may get trapped in long-lasting quasi-stationary states 
(QSS's), whose lifetime diverges with system size [3111]. 
Importantly, when performing the limit N — > oo (N be- 
ing the number of particles) , the system remains perma- 
nently confined in QSS's [5J |B]. As a consequence, for 
large long-range interacting systems, QSS's are directly 
accessible through experiment [7H3 ■ 

Until today, the large majority of studies aimed at elu- 
cidating the fundamental properties of LRIS's have been 
carried out on isolated systems, i.e. under the assump- 
tion that the system properties are not influenced by the 
external environment. However, recognizing whether a 
non-equilibrium QSS is stable to an external perturba- 
tion is of great importance [10] , both from a theoretical 
and an experimental point of view. A related funda- 
mental problem concerns the mechanism through which 
a LRIS exchanges energy with the surroundings. These 
questions epitomize the main motivation of the present 
work. 

The non-equilibrium dynamical properties of the 
LRIS's in contact with a thermal bath have been stud- 



ied for the first time only recently [TTT - 1T3] . As a possible 
realization of thermal bath, these authors considered a 
large Hamiltonian system with nearest-neighbor interac- 
tions, coupled to a fraction of the spins in the system. 
They concluded that the coupling with the bath intro- 
duces a new time scale in the evolution of the system: 
the weaker the coupling strength, the longer the system 
remains trapped in a QSS before relaxing to equilibrium. 

At variance with the above studies, we investigate here 
the dynamics of a LRIS in long-range contact with an ad- 
ditional large system trapped in a QSS. This interaction 
scheme can be regarded as a more clear-sighted realiza- 
tion of a thermal bath for a LRIS, opening the way to ap- 
plications in diverse fields such as cosmology and plasma 
physics. For example, one may think of the collision-less 
mixing between plasmas, or the operation of magnetic 
fusion devices for energy production or the merging of 
globular clusters to a self-gravitating galaxy. Further- 
more, it is also tempting to speculate that our simple 
scheme could be somehow relevant for the self-consistent 
interaction between dark (the bath) and baryonic (the 
system) matter in the universe (see, i.e. Ref. [2]). 

As a reference case, we have selected the Hamiltonian 
Mean Field (HMF) model [15], widely regarded as a pro- 
totype LRI system. The HMF Hamiltonian describes the 
one-dimensional motion of N rotators coupled through a 
mean field cosine-like interaction, 

N N 

where 0j is the orientation of the j-th rotator and pj its 
conjugated momentum. To monitor the evolution of the 
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system, it is customary to introduce the magnetization 
M, an order parameter defined as 



can be accomplished through the following normalization 
condition 
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where m, 
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The infinite-range coupling between rotators is respon- 
sible for the emergence of rather intriguing behaviors, 
including the existence of QSSs. In a QSS the system 
displays non-Gaussian velocity distributions and it takes 
values of M different than those predicted by equilibrium 
thermodynamics [HISIITt)]. 

Rigorous mathematical results |17j indicate that in the 
limit N — >• oo the discrete HMF dynamics reduces to the 
continuum Vlasov equation 
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where f(0,p,t) is the microscopic one-particle distribu- 
tion function, V(6)[f] = -M x [f]aos(6) - M y [f] sin(0), 
M x [f\ = Ho /(*,!>,*) cos 0d0 dp and M v\f\ = 
J-i T f™o O f( e >P> t ) sin0d9d P- The s P ecinc energy h[f] = 
J J(p 2 /2)f(9,p,t)d9dp - [Ml + Ml - 1) /2 is a con- 
served quantity. The Vlasov equation defines the nat- 
ural framework to address the puzzle of QSS's emer- 
gence [31 [5]. Specifically, QSS's are connected to the 
stable stationary solutions of the Vlasov equation. This 
observation suggests a statistical mechanics approach, in- 
spired by the seminal work of Lynden-Bell |18j , to char- 
acterize analytically the QSS properties. Lynden-Bell's 
approach is based on the definition of a locally-averaged 
( "coarse-grained" ) distribution, yielding an entropy func- 
tional defined from first-principle statistical-mechanics 
prescriptions. By constrained maximization of such an 
entropy, one obtains closed analytical expressions for the 
single-particle distribution in the QSS regime [5, 6 . As a 
natural consequence, the QSSs can be equally interpreted 
as equilibrium configurations of the corresponding con- 
tinuous description [19] . Hence, the QSS thermal bath 
that we consider here corresponds to a magnetized equi- 
librium solution of the underlying Vlasov equation ([3| . 

Let f B (9,p) be the normalized single-particle distribu- 
tion that characterizes the QSS bath. Such a function is 
obtained as the stationary solution of the Vlasov equa- 
tion ([3| corresponding to a water-bag initial distribu- 
tion, f (p,9) = l/[4A0 B A Ps ] for 9 € [-A9 B ,A9 B ],p E 
[— Ap B , Ap B ] and zero elsewhere. Note that the initial 
magnetization of the bath (Mq)b and its energy den- 
sity eb can be expressed in terms of A9b and Ap B , as 
(M Q ) B = smA9 B /A9 B and e B = Ap%/6 + 1/2(1 - 
(Mq) b ). This in turn implies that the initial water-bag 
profile is uniquely determined by (Mo) b and e B , in agree- 
ment with the Lynden-Bell theory [21] . 

At this point, t = in our discussion, another HMF 
system with water-bag profile is injected and let evolve 
consistently with the bath. This system, S in the follow- 
ing, is described in terms of its associated single-particle 
distribution fs(0,p). Clearly the bath should be signifi- 
cantly larger than the system to which it is coupled. This 



J f s (0,p,t)d0dp=l- J f B (9,p,t) d9dp 
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where e <C 1 sets the relative size of the two mutu- 
ally interacting S and B HMF systems. We are in- 
terested in tracking the time evolution of the distribu- 
tion f(9,p,t) = f B (9,p,t) + fs(9,p,t) under the con- 
straint Q. From the physical point of view, we are repro- 
ducing the microcanonical dynamics of one isolated HMF 
system (S + B), composed of two sub-systems supposed 
as distinguishable: the larger subsystem (the bath B) 
has already relaxed to its QSS equilibrium. The system 
S is initially confined in an out-of-equilibrium configura- 
tion of the water-bag type. To monitor the evolution 
of both subsystems, we follow the kinetic temperatures 
T a (t) = T a f p 2 f a (p,9,t)d9dp, with a = B,S and the 
corresponding magnetizations M a . Here, Ts = 1/e and 
Tb = 1/(1 — e). A typical time evolution of these observ- 
ables, obtained by numerical integration of the Vlasov 
equation ([3]), is illustrated in fig. [I] (see also supple- 
mentary material). After a short transient, the system 
reaches a quasi-equilibrium state where the mean value 
of the kinetic temperature is different from the temper- 
ature of the bath. In other words, the bath and the 
systems do not thermalize. Similarly, the two magne- 
tizations converge to different values. Importantly, we 
note that the specific values of temperature and mag- 
netization attained by the system spotlight a non-trivial 
interaction with the bath. Ts and Ms are indeed sub- 
stantially different from the values that the system would 
reach when evolved microcanonically from the same ini- 
tial condition. We obtain equivalent results upon sim- 
ulating the discrete A^-body dynamics ([!]). In this case, 
after a transient that gets progressively longer as the sys- 
tem size N = N s + N B is increased, AT = T B - T s 
and AM = M B — Ms tend to zero. Thus, granularity 
causes thermalization, which is instead prevented in the 
continuum (Vlasov) limit. We term canonical QSS's the 
quasi-equilibrium configurations that the system explores 
when in long-range contact with a QSS thermal bath in 
the zero energy-flux regime. 

In the continuum limit, when the system is trapped 
in a canonical QSS, we find that the average energy flux 
between the bath and the system indeed vanishes, mak- 
ing the two subsystems by all means decoupled and thus 
preventing thermalization (see supplementary material 
A for a more detailed analysis). It is remarkable that 
a zero-flux steady state is reached for T B ^ Ts in the 
non-collisional continuum limit, at variance with what is 
normally found in most collisional systems. 

Even more surprising is the behavior of the system dur- 
ing the "violent relaxation" stage towards the canonical 
QSS, which is characterized by a net energy flux from the 
(cold) bath to the (hot) system. To better illustrate this 
observation, we plot T B and Ts versus time in fig. [2} Note 
that Ts is larger than T B at t = 0, the time of injection. 
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FIG. 1: (Color online) Time evolution of temperature and 
magnetization. The bath QSS originates from a water-bag 
with energy 0.54 and initial magnetization 0.6. The system is 
initially homogeneous in space (zero magnetization) and its 
energy is set to 0.65. The coupling constant e = 0.024. 
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FIG. 2: (Color online) Bath and system temperatures ver- 
sus time. The system is initially space homogeneous (zero 
magnetization) and has energy 0.75. Inset: same plot with 
logarithmic time scale. Other parameters are as in fig. [I] 



As time progresses, the difference AT increases even fur- 
ther, resulting in an anomalous energy transfer from the 
bath to the system. In short, and counter-intuitively, 
the hot system gets hotter when placed in contact with 
a large long-range QSS reservoir. This observation, al- 
though fighting intuition, does not violate any laws of 
physics, as the second law of thermodynamics is only ex- 
pected to hold at thermal equilibrium. 

Once the system has settled down in its canonical QSS 
at zero average energy flux, AT and AM are found to be 
different from zero. In order to pinpoint the relation be- 
tween AT and AAf, we performed a series of simulations 
for the same bath conditions as specified in the caption 



of fig. [T] and varying the initial energy of the system S. 
Different energies lead to distinct canonical QSS's, as it 
happens to isolated systems trapped in microcanonical 
QSS's. At first glance, it is tempting to speculate that 
canonical QSS's might originate from a net balance of 
two opposing thermodynamic forces, presumably related 
to AT and AM. However, we find that the dynamical 
evolution of S is not influenced by the temperature of 
the bath Tb, at least for e <C 1, but only responds to 
its magnetization Mb- Therefore, provided Mb is kept 
fixed, Tb can be set to an arbitrary value, without signif- 
icantly altering the system dynamics. This is illustrated 
by the data collapse reported in fig. [3j 

This striking observation is unintuitive as compared to 
the case of short-range systems. Even more interestingly, 
it is by no means restricted to the HMF. In order to il- 
lustrate this fact, we note that in the Vlasov limit the 
distribution functions f a (9,p,t) (a — B,S) obey to 

-af+PW- v[fB + fs] W = Q (5) 

where V{9) is a generic mean-field potential (the prime 
defining ordinary differentiation with respect to 9), de- 
fined as 



V[f] = / f{9',p',t)v{9-9')d9'dp', 



(6) 



v(9 — 9') being the two-body potential. Since the sys- 
tem/bath relative size e <C 1, we can treat it as a per- 
turbative parameter, with fg ~ 0(e) and fg ~ 0(1). 
Expanding eqs. ^ and keeping only terms that cause 
changes ~ 0(e) in the physical observables, we are led to 
the two following coupled equations 
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The equation for the bath implies that this is frozen in its 
initial configuration, a stable equilibrium of the Vlasov 
equation, fg(q,p, t) = /s(<7,p> t = 0) at all times [22] , 
The equation for fg is the Liouville equation for a dis- 
tribution of uncoupled particles moving in an external 
potential, fs being constant. These conclusions are ut- 
terly general and should apply to any physical system 
whose density is governed by the Vlasov equation. For 
the HMF model, thanks to its rotational invariance, one 
has with no loss of generality 
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which is simply the Liouville equation for a set of un- 
coupled pendula. Hence, the leading-order evolution of 
fg depends only on Mb and not on Tb- Note that Mb 
sets the width of the resonance of the pendulum along p, 
which scales as ^/ Mb ■ This implies that the temperature 
should be proportional to Mg, as can be also appreciated 
by dividing eq. dsl) in the stationary state by y/Mg. 
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FIG. 3: (Color online) Difference between final and initial 
temperature of the system versus width of its initial water- 
bag, in reduced units. Data refer to different choices of the 
bath parameters and to different initial energies of the (ini- 
tially homogeneous) system. Symbols: direct integration of 
eqs. <[§ for Mb € [0.1,0.55], T B G [0.3,0.4]. Solid line: nu- 
merical solution of eq. B. Inset: Tsj vs. M$j for the same 
choice of parameters for the bath as in fig. [I] Circles: N- 
body simulations (average over 100 independent realizations), 
Nb = 4 x W :i ,Ns = 10 2 . Crosses: direct integration of the 
Vlasov equations. Solid line: integration of eq. nSD - 



Consistently with the above scaling arguments, we 
plot in fig. [3] (Tsj — Ts^)/Mb as a function of the 
rescaled width of the initial water-bag Ap$ / \J Mb, for 
different values of the bath magnetization and temper- 
ature. The data refer to direct integration of the (con- 
strained) Vlasov equations ([5| and of eqs. pi). In all 
cases, the data collapse nicely on a single master curve, 
which confirms the validity of our reasoning. An analyt- 
ical calculation of (Tsj — Ts,{)/Mb for Aps = yields 
(Tsj — Ts,i)/MB ~ 0.751, in excellent agreement with 
the result of direct integration of eq. ^ (see supplemen- 
tary material B and also Ref. [2U]). The inset further 
shows that iV-body simulations agree with all results ob- 
tained in the continuum limit. 

Summarizing, we have proposed an implementation of 
long-range QSS bath. We showed that a small system in 
true long-range contact with a large, long-range reservoir 
reaches a zero-flux steady state, that we term canoni- 
cal quasi-stationary state. These are stationary states 
of the system-bath coupled Vlasov equations, but quasi- 
stationary solutions of the associated V-body problem. 
Remarkably, in the explored range of parameters, we find 
that hotter-than-bath systems become hotter in canoni- 
cal QSS's. In the contest of the HMF model, based on 
simple scaling arguments, we have unveiled how the sys- 
tem anomalously increases its temperature as the frac- 
tion of its particles trapped in the resonance set by the 
bath magnetization gain energy. The kinetic energy gain 
is proportional to the value of Mb and independent of 



the bath temperature at the leading order in e. Fur- 
thermore, we argue that the results illustrated in this 
paper should be general and extend beyond the HMF 
case-study, whenever the collision-less Vlasov picture is 
a good description of the dynamics. Finally, our results 
might open the unexpected possibility of extracting en- 
ergy from a cold reservoir through long-range coupling. 

Dicussions during the workshop " Equilibrium and out- 
of-equilibrium properties of systems with long-range in- 
teractions" held at the Centre Blaise Pascal, ENS-Lyon 
in August 2012 are acknowledged. 



Supplementary Material A 

The energy flux from the bath B to the system S is de- 
fined as $b^s — —dEs/dt, where Eb is the total energy 
of B. In order to derive an explicit expression for 3>b^s, 
we start by calculating <f>j(t), the rate of energy loss of 
the j-the particle. Denoting by hj its energy, we have 



dhj(t) 
dt 



1 dM 

2~dt 



where nx, = (cos0,-,sin0j) and mj- = (— sin #7, cos 
M being the global magnetization 



1 N 

M = (M x ,M y ) = -J2 
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(9) 
(10) 



Here N is the total number of particles, i.e. the sum of 
those belonging to the bath, Nb, and those in the system, 
Ns- Summing over all particles belonging to the bath in 
eq. ([9]), one eventually obtains 

$B^ S = 22 ^ = -9 Eft M • m i + — — ■ Mb 



jeB 



2 dt 



(11) 

where Mb = J2jeB m j/^-B is the magnetization of the 
bath and the time derivative of M reads 



dM 

~dt 



1 N 



(12) 



3=1 



In the continuum limit the sums are replaced by integrals 

l - J P M- m ± (6)fB(p,8)dedp+^-M B 

(13) 



where 



M 



B 



f B m{0)d8dp M 



(f B +f s )m(6)d0dp 
(14) 

with m(9) = (cos 9, sin 0), m. ± (9) = (— sin 6*, cos 6*). 

According to the adopted sign convention, $B^s(t) is 
positive if the bath B cedes energy to the system. In 
fig. [4] the instantaneous energy flux (upper panel) is plot- 
ted versus time for a typical realization of the Vlasov 
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FIG. 4: (Color online) Time evolution of the instantaneous 
(top) and cumulated (bottom) bath-to-system energy flux. 
The system is initially space homogeneous and has energy 
0.75. Other parameters are as in fig. 1 in the paper. 



dynamics. After an initial transient, (pB^s(t) oscillates 
around zero, implying that the bath B and the system 
S have established a zero-average-flux dynamical equi- 
librium. This condition corresponds to the emergence of 
the canonical QSS. Furthermore, the net energy flux is 
positive, a fact that can be appreciated by looking at the 
evolution of the cumulated flux (see lower panel of fig.|4|. 
This implies a net transfer of energy from the bath to the 
system. 



Supplementary Material B 

The phase space of the pendulum is foliated in trajecto- 
ries with constant energy 



e = — — Mb cos ( 

nl/2 



(15) 



hence, p(9) = V% [e + Mb cos 9) L/ * . We want to discuss 
an analytic estimate of the quantity {Tgj — Tgi)/Mg for 
Aps — and for an initial homogeneous system, Mg(t = 
0) = 0. This calculation has the merit of enabling one 
to gain insight into the nature of the canonical QSS and 
further clarify the scaling adopted in fig. 4. This analysis 
can be extended to cover the case Aps ^ 0, and also 
Ms(t = 0) 7^ 0, a generalization to which we shall return 
in a separate contribution. 

We note that T Sil = for Ap s = 0. To evaluate T SJ , 
we first consider the average kinetic temperature of the 
particles which are assigned a given energy e. In formulae 



(P 2 )e = 



T(e) 
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e 2 dt 
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where (-) e indicates a time average over one period 
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1Mb V 2M b 
K(-) being the complete elliptic integral of the first kind. 

Expression ( 16 1 takes the equivalent form 
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where 9(e) = cos _1 (— e/Ms) is the angle of inversion 
of the selected (closed) trajectory. By performing the 
integral one eventually gets 
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where E (•,•) is the incomplete elliptic integral of the sec- 
ond kind. The final temperature of the system can now 
be evaluated as 



m e 



t sj = (P ) = / (P ) e p(e)de 

J -Mb 



(20) 



where p(e) is the density of states of the system, which 
is univocally fixed by the initial condition. The integral 
in eq. (20) extends from —Mb to Mb, i.e. the energies 



that identify the separatrix of the pendulum. In fact, the 
system is trapped inside the separatrix |e| = Mb, given 
the specific condition selected here (Ap$ — 0, hence no 
particle lies outside the resonance at t = 0). Recalling 
eq. (151, the distribution p(e) can be calculated easily, as 



P(e) 



1 1 



(21) 



eventually obtains 



Plugging eq. (21) into eq.(20) and recalling eq. (19), one 
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Numerical integration gives T$j/Mb ~ 0.751, in excel- 
lent agreement with the data reported in fig. 3. In the 
general case (Aps ^ 0), e cx Ap s . The scaling suggested 
by eq. (22) implies Apl \[Mb, which in turn explains the 
origin of the reduced variables used in fig. 3. 
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